GLOBAL DISSIPATIVE DYNAMICS OF 
THE EXTENDED BRUSSELATOR SYSTEM 
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Abstract. The existence of a global attractor for the solution semifiow of the extended Brussela- 
tor system in the phase space is proved, which is a cubic-autocatalytic and partially reversible 
reaction-diffusion system with linear coupling between two compartments. The method of group- 
ing and re-scaling estimation is developed to deal with the challenge in proving the absorbing 
property and the asymptotic compactness of this typical multi-component reaction-diffusion sys- 
tems. It is also proved that the global attractor is an {H, E) global attractor with the L°° 
regularity and that the HausdorfE dimension and the fractal dimension of the global attractor 
are finite. The results and methodology can find many applications and further extensions in 
complex biological and biochemical dynamical systems. 



1. Introduction 

Consider an extended Brusselator system consisting of cubic-autocatalytic and partially re 
versible reaction-diffusion equations with linear coupling between two compartments, 
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diAu + a- {b + k)u + u^v + Di{w - u) + Nip, (1.1) 

d2Av + bu- u^v + D2iz - v), (1.2) 

d3Aip + ku-{X + N)ip + D3{i^-p), (1.3) 

diAw + a- {b + k)w + w^z + Di{u - w) + Nip, (1.4) 

d2Az + bw-w^z + D2{v- z), (1.5) 

d3AiJj + kw-{X + N)',jj + D3{ip-TP), (1.6) 



for {t,x) e (0,oo) X n, where fl C 5R"(n < 3) is a bounded domain with a locally Lipschitz 
continuous boundary, and the coefficients d[s,D[s,{i — 1, 2, 3), a, 6, fc. A, and N are positive con- 
stants. Note that there is a linear coupling between the two subsystems of components (m, w, ip) and 
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(w,z,ip), which can be interpreted as the interaction between two cells or compartments. Given 
the homogeneous Dirichlet boundary condition 

x) = v{t, x) — (p{t, x) — w{t, x) — z{t, x) — ipit, a;) = 0, i > 0, x G dVt, (1-7) 

and an initial condition 

u(0, x) = uo{x), u(0, a:) = vo{x), ip{Q, x) = ipo[x), 

it;(0, a;) = wq{x), z(0, x) — 2:0(2;), "0(0, ^) — i^o{x), x d 

we shall study the asymptotic dynamics of the semiflow of the weak solutions to this initial-boundary 
value problem pTT]) - (fTSl) . 

Historically the Brussels school led by the renowned physical chemist and Nobel Prize laure- 
ate (1977), Ilya Prigogine, made remarkable contributions in the research of complexity of cubic 
autocatalytic reactions. The mathematical model signifying their seminal work was named as 
the Brusselator, which is originally a system of two ordinary differential equations shown below, 
cf . [H[5B1I3S] . The scheme of chemical reactions described by the Brusselator kinetics is 

A — >\J, B + U — ^V + D, 
2U-I-V — >3U, U — ^E, ^^'^^ 

where A, B, U, V, D, and E are chemical reactants or products. Let u(t,x) and v{t,x) be the 
concentrations of U and V. Assume that the concentrations of the input compounds A and B are 
held constant during the reaction process, denoted by a and b respectively. Also assume that these 
elementary reactions are not reversible. Then by the Pick's law and the law of mass action, and 
through non-dimensionalization, one can formulate a system of two nonlinear reaction-diffusion 
equations called (diffusive) Brusselator equations, 

— = d\l\u + a — (0 -I- k)u V, 

f (1.10) 

aw , . , 2 

— = d2i\v + bu — u V. 

at 

The known examples of autocatalytic reactions which can be modeled by the Brusselator equa- 
tions include ferrocyanide-iodate-sulphite reaction, chlorite-iodide-malonic acid (CIMA) reaction, 
arsenite-iodate reaction, many enzyme catalytic reactions, cf. [Tl[51|5]. 

In 1993 two consecutive papers [22^, 24 discovered a variety of interesting self-replicating pattern 
formation associated with cubic-autocatalytic reaction-diffusion systems, respectively by an experi- 
mental approach and a numerical simulation approach. Since then numerous studies by mathemat- 
ical and computational analysis have shown that the cubic-autocatalytic reaction-diffusion systems 
such as Brusselator equations and Gray-Scott equations [14lil5j exhibit rich spatial patterns (includ- 
ing but not restricted to Turing patterns) and complex bifurcations [T]|5ll6ll9lfT0lll2lll6lll7ll20ll25ll27ll38j 
as well as interesting dynamics [8l fTTl[T8l[2n[28H30l[39ll40] on ID or 2D domains. 

For Brusselator equations and the other cubic-autocatalytic model equations of space dimension 
rt < 3, however, we have not seen substantial advancing results in the front of global dynamics until 
recently [iTHiS] . 

The extended Brusselator system in this work features the coupling of two compartment variables 
and the reversible reaction U E in both compartments, where E is the main product being 
removed at a constant rate A. This system consisting of six coupled components with partial 
reversibility is closer to the biological network dynamics |13lfT^ . 
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Multicell models in biological dynamics generically mean the coupled ODEs or PDEs with large 
number of unknowns representing component substances, which appear widely in the literature of 
systems biology as well as cell and molecular biology [32 1 133 1 137 ] . 

For this extended Brusselator system and many other biochemical activator-inhibitor reactions, 
the asymptotically dissipative condition in vector version 

lim /(s) ■s<C, 

where C is a nonnegative constant, is inherently not satisfied by the nonlinearity involving the 
opposite-signed autocatalytic terms and the opposite-signed coupling terms. This and the addi- 
tional equations of variables ip and ip due to the partially reversing reactions cause substantial 
difhculties in a priori estimation of all six components in proving the absorbing property and the 
asymptotic compactness of the solution semiflow. The novel mathematical feature in this paper 
is to overcome these obstacles by a method of grouping and re-scaling estimation. This method- 
ology will be potentially very useful in analyzing biological network dynamics of multicomponent 
reaction-diffusion systems. 

We start with the formulation of an evolutionary equation associated with this system. Define 
the product Hilbert spaces as follows. 

These six-component spaces are the phase spaces of different regularity for the component functions 
u{t, ■),v{t, ip{t, ■),w{t, ■), z{t, ■) and ip{t,-). The norm and inner-product of or any component 
space will be denoted by || • || and (•,•), respectively. The norm of L^'(f2) or the product 

space L''(r2) = [L''(f2)]^ will be denoted by || • ||lp if p 7^ 2. By the Poincare inequality and the 
homogeneous Dirichlet boundary condition (|1.7|) . there is a constant 7 > such that 

||Vef>7l|^f, forUH^movE, (1.11) 

and we take ||V^|| to be the equivalent norm ||^1|_e of the space E or the equivalent norm ||CII_ffi(n)- 
We use I • I to denote an absolute value or a vector norm in a Euclidean space. 

By the Lumer-Phillips theorem and the analytic semigroup generation theorem [34 , the linear 
operator 

MA 
dzA 
^ dsA 
~ diA 

V 

is the generator of an analytic Cg-semigroup on the Hilbert space H, which will be denoted by 
{e'**}t>o. By Sobolev embedding theorem, Hq^^I) ^ L^(r2) is a continuous embedding for n < 3. 
Using the generalized Holder inequality, we have 

\\uM < hllielklUe, \\W^Z\\ < \\w\\le\\z\\Le, U,V,W, Z ^ L^Q). 






\ 




















d2A 








dsA/ 



■.D{A){=n) — >H (1.12) 
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E 



(1.13) 



By these facts one can verify that the nonlinear mapping 

f a~{h + k)u + u^v + Di{w - u) + N(f\ 
hu — v?v + D2{z ~ v) 
ku- {\ + N)lp + D:i{i} - tp) 
- (6 + k)w + w'^z + Di{u - w) + Nijj 
bw — VlP'Z + D2{v — z) 
\ kw - {X + N)i' + Dsi^ - i;) ) 

where g = {u,v,ip,w^ z^ip), is a locahy Lipschitz continuous mapping defined on E. Then the 
initial-boundary value problem (|l.ip - (|1.8p of this extended Brusselator system is formulated into 
an initial value problem: 



^^Ag + fig), t>0. 
5(0) = 90 = col{uo,VQ,ipQ,wo,zo,ipo) e H. 



(1.14) 



Here g{t) = col {u{t),v{t),(p{t),w{t),z{t),il;{t)) can be written as {u{t),v{t),ip{t),w{t), z{t),^jj{t)) 
for simplicity. Accordingly we shall write go = {uQ,vo,(po,wo, zo,ipo)- 

The following proposition will be used in proving the existence of a weak solution of the initial 
value problem (|1.14p . Its proof is seen in [71 Theorem II. 1.4] and in [4j Proposition 1.3.3]. 



Proposition 1.1. Consider the Banach space 

WiO, r) = {tj{-) : 77 G L^{0, r; E) and dtij G L^iO, r; E*)} 
with the norm 

WvWw = \\v\\l^O,t-E) + \\dtV\\L^O,T-E-)- 

Then the following statements hold: 

(a) The embedding W{0,t) ^ L^(0,t; H) is compact. 

(b) If rj G W(0,t), then it coincides with a function in C([0,t];7J) for a.e. t G [0,t]. 

(c) If rj,C, G W{Q^t), then the function t {rj{t) ^ C^{t)) h is absolutely continuous on [0,t] 



(1.15) 



and 



d 
dt 



drj 
'dt' 



dc 

dt 



(a.e.)tG [0,t], 



where (•, •) is the {E*,E) dual product. 



By conducting a priori estimates on the Galerkin approximations of the initial value problem 
(|1.14p and the weak/weak* convergence, we shall be able to prove the local and then global existence 
and uniqueness of the weak solution g(t) of (|1.14p in the next section, also the continuous dependence 
of the solutions on the initial data and the regularity properties satisfied by the weak solution. 
Therefore, the weak solutions for all initial data in H form a semiflow in the space H. 

We refer to [THHUfM] and many references therein for the concepts and basic facts in the theory 
of infinite dimensional dynamical systems, including few provided here for clarity. 

Definition 1.2. Let {S{t)}t>o be a semiflow on a Banach space X. A bounded subset Bq of X 
is called an absorbing set in X if, for any bounded subset B C X, there is some finite time tQ > 
depending on B such that S{t)B C Bq for all t > to. 
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Definition 1.3. A scmiflow {S{t)}t>o on a Banach space X is called asymptotically compact in X 
if for any bounded sequences {xn} in X and {t„} C (0, oo) with t„ — > oo, there exist subsequences 
{xrik} of {xn} and of {t„}, such that linifc^oo S{tn^)xnk exists in X. 

Definition 1.4. Let {S{t)}t>o be a semiflow on a Banach space X . A subset of X is called a 
global attractor in X for this semiflow, if the following conditions are satisfied: 

(i) ^ is a nonempty, compact, and invariant set in the sense that 

S{t)£/ = £/ for any t > 0. 

(ii) £/ attracts any bounded set B of X in terms of the Hausdorff distance, i.e. 

dist(S'(t)J5, £/) = sup inf \\S{t)x - y\\x ^0, as i oo. 

Here is the main result of this paper. We emphasize that this result is established unconditionally, 
neither assuming initial data or solutions are nonnegative, nor imposing any restriction on any 
positive parameters involved in the equations ()1.1|) - ()1.6|) . 

Tiieorem 1.5 (Main Theorem). For any positive parameters di,d2,d3, Di, D2, D3,a,b,k, X, and 
N , there exists a global attractor si in the phase space H for the semiflow {S{t)}t>Q of the weak 
solutions generated by the extended Brusselator evolutionary equation (|1.14p . 

The following proposition states concisely the basic result on the existence of a global attractor 
for a semiflow, cf. [31,34,36 . 

Proposition 1.6. Let {S(t)}t>o be a semiflow on a Banach space X. If the following conditions 
are satisfied: 

(i) {S'(i)}t>o has an absorbing set in X , and 

(ii) {S{t)}t>o is asymptotically compact in X , 

then there exists a global attractor s/ in X for this semiflow, which is given by 

s/ = i,iBo)= fl Clx \J{S{t)Bo). 

T>0 t>T 

In Section 2 we shall prove the global existence of the weak solutions of the extended Brusselator 
evolutionary equation (jl.l4l) . In Section 3 we show the absorbing property of this extended Brus- 
selator semiflow. In Section 4 we shall prove the asymptotic compactness of this solution semiflow. 
In Section 5 we show the existence of a global attractor for this semiflow and its properties as being 
the (7J, E) global attractor and the L°° regularity. We also prove that the global attractor has 
finite Hausdorff dimension and fractal dimension. 

2. Global Existence of Weak Solutions 

We shall write ult, x),v{t, x), (p(t, x), w(t, x), z{t, x) and ^p{t, x) simply as u{t),v{t), Lp{t), w{t), z{t) 
and '4'{t), or even as u, v, tp^ w, z and ip. Similarly for other functions of {t, x). 

The local existence of the solution to a system of multicomponent reaction-diffusion equations, 
such as the initial value problem (IVP) of (|1.14p . with certain regularity requirement is not a trivial 
issue. There are two different approaches to get a solution. One approach is the mild solutions 
provided by the "variation of constant formula" in terms of the associated linear semigroup {e^*}t>o 
but the parabolic theory of mild solutions requires that go ^ E instead oi go E H assumed here. 
The other approach is the weak solutions obtained through the Galerkin approximations and the 
Lions-Magenes type of compactness treatment, cf. [71 [^[5Tj . 
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Definition 2.1. A function g{t,x), {t,x) G [0,t] x ft, is called a weak solution to the IVP of the 
parabolic evolutionary equation (|1.14l) , if the following two conditions are satisfied: 

(i) -^{g, r]) = {Ag, 77) + (/(g), 77) is satisfied for a.e. t G [0, r] and any i] ^ E; 

(ii) g{t, •) e L'^{0,t;E) n Cu,([0, r]; iJ) such that 5(0) = go- 
Here (■, ■) stands for the (E* , E) dual product. 

For reaction-diffusion systems with more general polynomial nonlinearity of higher degrees, the 
definition of corresponding weak solutions is made in |71 Definition XV. 3.1]. 

Lemma 2.2. For any given initial data go £ H, there exists a unique local weak solution g{t) = 
{u{t),v{t),(p{t),w{t), z{t),ip{t)), t 6 [0,t] for some t > 0, of the IVP of the extended Brusselator 
evolutionary equation (|1.14p . which satisfies 

g e C([0, T];H)n C\{0, r); ) n L^O, r; E). (2.1) 

Proof. Using the orthonormal basis of eigenfunctions {ej{x)}^^ of the Laplace operator with the 
homogeneous Dirichlet boundary condition: 

Acj + XjCj =0 in fl, Cj \gn= 0, j = 1, 2, • ■ • , 

we consider the solution 

m 

= 5I«rW^j(^)' t^[0,r],xen, (2.2) 

of the approximate system 

%^ = Agm + Pmfigm), t>0, 

at (2.3) 

5m (0) = Pm.90 e 

where each q^{t) for j = 1, • • • , m is a 6-dimensional vector function of t corresponding to the six 
unknowns u, v, yj, w, z and ■0, and Pm ■ H Hm = Spanjei, • • • , Cm} is the orthogonal projection. 
Note that for each given integer to > 1, (|2.3p can be written as an IVP of a system of ODEs, whose 
unknown is a Gm-dimensional vector function of t. The components of the solution to this system 
of ODE are the coefficient functions in the expansion (|2.2I) oi g,n{t,x), namely, 

g™(t) = col (97jt), g™(t), g7;(i), g™ (t), g;-':(t), J = 1, • • • , m). 

The IVP of this system of ODE can be written as 



^ = A„g™(i) + F™((7"(i)), t > 0, 
at 



(2.4) 



5^(0) = Col(P™Uoj, PmVOj, Pm^Qj, PmWQj , PmZQj, PmlpOfi j = 1, " ' ' , m) . 

Note that is a matrix and i^™ is a 6TO-dimensional vector of cubic-degree polynomials of 6m- 
variables, which is a locally Lipschitz continuous vector function in M^™. Thus the solution of the 
system of ODE ()2.4p exists uniquely on a time interval [0,t], for some r > 0. Substituting all the 
components of this solution q"^{t) into (|2.2p . we obtain the unique local solution gm{t, x) of the IVP 
TO > 1. 

By the multiplier method we can conduct a priori estimates based on 

^llffm(i)llL + (dV5™(0, Vg™(0)if„ = / (5™ (i)), 5™ (*))//„, t e [0,t], 
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where the matrix d — diag((ii, (^2, da, di, (i2, (is). These estimates are similar to what we shah 
present in the next lemma. Note that ||(?m(0)|| = ||i^m.9o|| < ||.9o|| for all m > 1. It follows that (see 
the proof of the next lemma) 

{5ni}m=iis a bounded sequence inL^(0,T;£;) f) L°°{0,t; H), 

so that 

{^5m}m=iis a bounded sequence mL^{0,T; E*), where E* = [H^'^ (fl)]'' . 

Since f : E ^ H is continuous, 

{Pnif{gm)}^=iis a bounded sequence in 2.^(0, r;i?) C L^{0,t;E*). (2.5) 

Therefore, by taking subsequences (which we will always relabel as the same as the original se- 
quence), there exist limit functions 

g{t,-) e L^iO,T;E)nL°°{0,T;H) and $(t, •) e L^(0, r; i?) (2.6) 

such that 

gm — > g weakly in £^(0, t; E), as to — oo, 

gm — > g weak* in L°°(0, r; H), as m oo, (2.7) 
^ffm — ^ -^g weakly in i^(0, r; i?*), as to — oo, 

and 

PrnfiSm) — > weakly in L^(0, r; H), as to -j- oo. (2.8) 

To estimate the distributional time derivative sequence {dtgm}m=ii take the {E*,E) dual prod- 
uct of the equation (|2.3p with any 77 £ £' to obtain 

\\dtgrn{t)\\E' < \\Agrn{t)\\E' +c\\Pmf{gAt))\\H, t € [0,t], TO > 1, 

where c is a universal constant, so that by the aforementioned boundedness it holds that 

{dtgm}m=i^^ ^ bounded sequence inL^(0, r; i?*). 

From extracting a convergent subsequence and by the uniqueness of distributional derivative it 
follows that 

dtgm — > dtg weakly in L^{0, r; E*), as m ^ oo. (2.9) 
In order to show that the limit function 5 is a weak solution to the IVP of (|1.14p . we need to 
show $ — f{g). By item (a) of Proposition [TTTl the boundedness of {gm} in L^{0, t; E) and {dtgm} 
in i^(0, r; E*) implies that (in the sense of subsequence extraction) 

g,n — > g strongly in £^(0, r; iJ), as ni 00. (2-10) 
Consequently, there exists a subsequence such that 

gm{t, x) — > g{t, x) for a.e. {t, x) £ [0, r] x fi, as m — )• 00. (2-11) 
Due to the continuity we have 

Pmfigmit, x)) — > f{g{t, x)) for a.e. {t, x) G [0, r] x fi, as to -J> 00. (2.12) 
According to 23, Lemma 1.1.3] or [7, Lemma II. 1.2], the two facts (|2.5p and (|2.12[) guarantee that 
Pmfigrn) — > f {g) wcakly in L'^{0,t;H), as to 00. (2.13) 
By the uniqueness, (|2?8l) and ([2J3l) imply that $ = f{g) in L^(0,t; H). 
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With (|2.7p . (|2.9p and (|2.13p . by taking limit of the integral of the weak version of (|2.3p with any 
given 77 G i^(0, r; E), 

■ ri) dt = [{Ag,n,ri) + {Pmfigm), v)] dt, as to-> 00, 



we obtain 

'dg 



V dt 

V]dt^l [{Ag,rj) + {f{g),rj)]dt, for any 7jGL-'{0,T- E). (2.14) 



dt 

Let 77 S be any constant function. Since in (|2.14p one can replace [0, r] by any small time interval, 
it follows that 

^(9,V) = {A9,V) + ifi9),v), for a.e. t e [0,t], and for any 77 e (2.15) 

Next, dtg G L'^{0,t;E*) shown in (gH) implies that g e C^{[0,t]; E*). Since the embedding 
H ^ E* is continuous, this g G Cw{[0, t];E*) and g G L°°(0, r; i?) shown in (|2.7p imply that 

gGC^([0,T];iJ), 

due to f?! Theorem II. 1.7 and Remark II. 1.2]. Now we show that ^(O) — go. By item (c) of 
Proposition ll.il for any 77 G C^{[0,t]; E) with 77(t) = 0, we have 

i-g, dtv) dt^ f [{Ag, 7;) + (/(g), 7;)] dt + (g(0), 77(0)), 
Jo Jo 

and 

/■T fT 

{-9m,dtri) dt= [{Ag^,T]) + {P„if{g„i),T])]dt+{P„igo,v{^))- 
Jo 

Take the limit of the last equality as 7n — > 00. Since Pmgo ^ 9o in H, we obtain (g(0),77(0)) = 
{go,'il{0)) for any 77(0) G E. Finally, by the denseness of E in H, it implies that g{0) — go in H. 
Then by checking against Definition 4, we conclude that the limit function g is a weak solution to 

the ivp of (mm). 

The uniqueness of weak solution can be shown by estimating the difference of any two possible 
weak solutions with the same initial value go through the variation of constant formula and the 
Gronwall inequality. 

By item (b) of Proposition ll.il and the fact that the weak solution g G W{0, r), the space defined 
in (jl.lSp . we see that g G C{[0,t]; H), which also infers the continuous dependence of the weak 
solution on go. 

Moreover, since g G L'^{0,t; E), for any t G (0,t) there exists an earlier time to G {0,t) such 
that g{to) G E. Then the weak solution coincides with the strong solution expressed by the mild 
solution on [io5''"]i cf. [7ll34| . which turns out to be continuously difFerentiable in time with respect 
to iJ-norm for t G (^O: ■''): cf. [34j Theorem 48.5]. Thus we have shown g G C^((0, r); H). Therefore, 
the weak solution g satisfies all the properties in (|2.ip . The proof is completed. □ 

Lemma 2.3. For any initial data go — {uq, vo, fo, Wq, zq, ipo) £ H , there exists a unique global weak 
solution g{t) = {u{t)^v{t),ip{t),w{t),z{t),il){t)), t G [0, 00), of the IVP of the extended Brusselator 
evolutionary equation (|1.14p and it becomes a strong solution on the time interval (0,oo). 
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Proof. Taking the inner products ( (|1.2I) . v(t)) and ( (jl.Sp . z(t)) and then summing them up, we get 

= / {-u^v^ +buv -w'^z'^ + bwz - D2[v'^ -2vz + z'^]) dx (2.16) 

= [ - {{uv - h/2f + [wz ~ b/2f + D2{v - zf) dx + ]-b^\n\ < ^h^^. 
Jn 2 2 

By Poincare inequahty, it follows that 

j^{\\vr+\\zr)+2jd2{\\vr+\\zr)<b'\ni 



which yields 



\Ht)f + \\z{t)f < e-^-""-' (||«of + \\zof) + for t G [0,T„,.). (2.17) 

27^2 

In order to treat the (m, i^)-components and the (w, '(/')-components, we shall go through several 
steps. Let yit, x) = u{t, x) + v(t, x) + w(t, x) + z{t, x) and ^(t, x) = ip{t, x) + ipit, x). Add up the 
equations (|l.ip . (|1.2I) . (11.41) and (II. 5p and add up (|1.3I) and (II. 6p . respectively, to get the following 
equations satisfied by y{t, x) and ^{t, x), 

dv 

= diAy -ky+ [(rfa - di)A{v + z) + k{v + z) + 2a] + N^, 

^ = dsA^ + ky~ k{v + z)~{\ + N% 
Rescaling ^ — jiS with /i — k/N, we get 

^ = diAy -ky+ [(ds - di)A{v + z) + k{v + z) + 2a] + kS, (2.18) 

^— = ^idsAS + ky- k{v + z) - (^A + k)~. (2.19) 
Now taking the inner-product ( (|2.18p . vit)) we obtain 

\i:\\vV + rfillVyf < / [(d2 - d,)A{v + z) + {v + z) + 2a]ydx - fc||yf + k\\v\\\\S\\ 



2dV 

< Ml - rf2|||V(« + z)||||Vy|| + \\v + z\\\\y\\ + 2a\nny\\ ~ fc||yf + k\\y\\\\S\\ 

^ ^ 2 ^ + ,)||2 + ^|,y||2 + 2 8 

4 di 4 di-f dij 



k\\yr + k\\y\m 

,MlVyf + ^ ^ , . , - 

2 di ai7 ai7 



< ^llVyf + Ml^||V(. + z)f + + + 3^a^|l^| - kWyf + fc||y|||| -||, 



so that 



ihf + A II v„f < ^^11 V(. + .)||^ + (||..f H- ||.f + 2.'m) ^^^^^ 

-2t||af + 2t||!,||||E||, 
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Then by substituting (PTf]) into (E?^ wc get 

j^\\y\? + rfillVyiP < ^^^^^l|V(^; + z)f - 2k\\yf + 2A;||y|| ||H|| + Ci(z;o,zo,0, (2-21) 
where 

16a2 



o / 41,2 

Ci(z;o,zo,t) = ^e-^^^^^* (||i;of + ||zof ) + [7^^ + 



Taking the inner- product ( (|2.19p . S'(i)) we obtain 



(2.22) 



/.-ll^f + 2Aid3|| V^IP < 2fc||i; + z|| ||S|| + 2fc||y|| ll^ll - 2(/.A + fc)||Hf 

< A^d37ll- f + + 2fc||y|| ||H|| - 2(mA + fc)|| -IP, 

where ^(i37||S'|p < /xd3||VH|p, so that 

+ MdsllVHf < 2fc||y|||| -II - 2(mA + fc)|| + C^ivo^z^^t), (2.23) 



in which, by (PTT)) . 

C2(W0, 2^0,0 



^^37 



e -'-'^ u;o r+ ^0 



(2.24) 



(2.25) 



27^2 

Then adding up p.2ip and ()2.23p and noticing that 

-2fc||yf + 4fc||j;|||| -II - 2(mA + fc)|| -f = -2 [fcdiyll - 1| -[D^ + ;.A1| -f ] < 0, 

we obtain 

<^^i^||V(. + z)f + 

where Ci and C2 are given by (|2.22l) and (|2.24p . Integration of the inequahty (|2.25p with the term 
niin{(ii, (i3}(|| Vyjp + IIa'~^^'^VCIP) dropped yields the foUowing estimate, 

\\y[tW + WmV < ""^""ll'^Iil (11^0 + + ^0 + ^olP + 11^0 + M') 
min{l,/i ^} 

^^^1^^ r||V(.(.) + z(.))fd. (2-26) 



min{l,/z i}di Jo 

+ C3 {\\vof + \\zof) + C4\fl\t, te [0,T,nacc), 

where we get from (|2.22l) and (I2.24p that 

3— + —r- ' = niax{l, Ai} 



C3 



27^2 V'^17 t^dal J ' ' V7^'^ic^2 rfi7 fJ-l^d2d3 

From dmi]), for t G [0,r™aa;) we have 

d2 f \\V{vis)±zis))\\'ds<2d2 r(||Vi;(s)f + ||Vz(s)f) ds < (||^;of + H^of ) +6'|^^|^. (2.27) 
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Substitute this into (|2.26p to obtain 

\\y{t)r + wmr < (ii^o +v„+wo+z,\\' + yo+ m') 

nnn|l,/i ^| 

+ ( -^It'-uL + ^0 (11"°"' + 11^°"') (2.28) 
\min|l,^ |ai«2 / 



2'^^ ' C,]\n\t, te[0,T_.). 



^min{l,/i ^}(iifi2 
From (PT7|) and (P?^ it follows that 

\\u{t) + w{t)r + Mt) + mr = im («w + + ii^w + mr 

< 3(||2/(t)f + \\vit)r + P(<)f ) + Ur < CsdlffolP + t e [o,r„_), 



(2.29) 



where C5 > is a constant independent of go and t. 

On the other hand, let p{t, x) = u(t, x) + v(t^ x) — w{t, x) — z(t, x) and 0{t, x) = {pit, x) — ip(t, x), 
which satisfy the equations 

^ = diAp -ik + 2Di)p + {d2 - di)A{v -z) + {k + 2(Di - D2)){v - z) + NO, 

— = d^M + kp- k{v - 2) - (A + N))e - 2D39. 
at 

Rescaling 6 = ^0 with the same /i as above, we get 

^ = diAp -{k + 2Di)p + {d2 - di)Aiv -z) + {k + 2(Di - D2))iv - z) + kO, (2.30) 
do 

H— = ^idsAO + kp- k{v - z) - (AiA + k))e - 2^iD30. (2.31) 
Taking the inner-product ( (|2.30p . v(t)) we obtain 

~\\pr+d^\\vpr+{k+2D,)\\pr 

< [di - d2)\\\7{v - z)||||Vp|| + \k + 2{D, - D2)\\\v - z\\\\p\\ + 

■ ^ \k+2{D,-D2m\v-zr+k\\e\\\\pi 



AD 



so that 



4lbf + rfillVpf + D,\\pf < "'^l' \\W{v - z)f + C6(^o, ^0,0 , , 

fli (2.33) 



2fciH|^ + 2fc||eiiiH|, 



where 



Ceivo,zo,t)^^\k + 2iDi-D2)\^ 



e-'-'^H\\vor+\\zor) + ^^^ 



(2.34) 
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Taking the inner-product ( (|2.3ip . 0{t)) we get 

M^lief + Aif^sllVef < 2k\\v - z\\\\0\\ + 2fc||p||||0|| - 2{nX + k + 2^Z?3)||6)f 

< ^'^""f n"'"'^ + 2^lbll He'll - 2(mA + k)\\ef, 

so that 



fi[j^\\0r + d3\\V0r + X\\er] <2k\\p\\\\O\\-2k\\0r + Cr{vo,zo,t), (2.35) 



where 

k _9-,,^„f/ii ii9 II ii9^ b 



C7{vo,Zo,t) 



g-27rf2tr||„„||2 j_ ||^„||2^, 



27^2 



(2.36) 



(2.37) 



2fiD3 

Summing up p.33p and (j2.35p and noticing that 

-2fc|b||2 + 4fc||e||||p|! - 2fc||^|l2 = -2fc(|b|| - ||©||)2 < 0, 

we have 

|(l|pf + h^-'Or) + min{di,d3}(||Vpf + U-'Vef) 
< ^^-7^I|V(« - z)\\' - {D.WpW' + M-^A||0f ) + Ce + Cj, 

where Cg and Cy are given by (ESI and (I23S1). We drop the term min{di, dgKH Vpf + ||/^-iV0f ) 
from the left side and the term — (_Di||pp + /x^^A||0p) from the right side of (|2.37p . Then integrate 
the resuhing inequahty and use (|2.27l) to get the fohowing estimate, 

niin|l,/i 



(2.38) 



+ niax{l,^}|di-d2r / ||V(w(s) -z(s))l|'ds 
«i Jo 

< "^^^^^'^ I (11^0 + _ ^ + ||(/?o - V-of ) 
min|l,/i ^1 



where 



max{l,/x}|di - d2|- , max{l,^} fc^ 1 n 

c(ic(2 27d2 \2iJ,D3 Di ) 

From (PTT)) and ((051) it follows that 

\\u{t) - w{t)f + ii^(t) - mf = \w) - (vit) zit))r + Mt) - mr 

<3 (Ib(t)f + \\v{t)\\' + \\z{t)\\') + \\e{t)r < CgdLgof + \m t e [0,T„,„,), 
where Cg > is a constant independent of initial data go and time t. 



(2.39) 
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Finally, from (I2.17p . (|2.29p and (|2.39p we conclude that for any initial data go G H, all the six 
components of the weak solution g{t) — g{t]go), i.e. v{t), z(t), and 

- knt) + w{t)) + Ht) - wim, wit) = huit) + wit)) - iuit) - wit))], 

f i (2.40) 

m - 2 w + ^ w) + w - V' w)] , m = 2 iivit) +m)- (At) - m)] , 

are uniformly bounded in [O^Tmax) if Tj^ax is finite. Therefore, for any S H, the weak solution 
git) of (|1.14p will never blow up in H at any finite time. It shows that for any G H, Tmax = oo, 
and the weak solution exists globally. Moreover, as shown in the last paragraph of the proof of 
Lemma 12.21 any weak solution turns out to be a strong solution on the time interval (0, oo). The 
proof is completed. □ 



3. Absorbing Properties 



By the global existence and uniqueness of the weak solutions and their continuous dependence 
on initial data shown in Lemma 12.21 and Lemma 12. 3[ the family of all the global weak solutions 
{git; go) ■ t > 0,go G H} defines a semifiow on H, 

Sit) : go git; go), go^H, t>0, 

which will be called the extended Brusselator semifiow associated with the extended Brusselator 
evolutionary equation (jl.l4l) . 



Lemma 3.1. There exists a constant Ki > 0, such that the set 

Bo^{g^H:\\gf<K^} 
is an absorbing set in H for the extended Brusselator semifiow {Sit)}t>o- 
Proof. For the components and z(f) of this semifiow {5(t)}f>o, from (|2.17p we obtain 

limsup(||«(Of + \\zit)f) <Ro = 



(3.1) 



(3.2) 



and for any given bounded set i? C iJ there exists a finite time T]^'^ > Osuchthat ||i'(t)|p + ||z(t)|p < 
Ro for any go G B and all t > Tg" . For any t > 0, (|2.16|) implies that 



'''\\\vvis)r + iivz(.)f ) ds < U\Ht)r + imr + b'm) 

£12 



< 



'27^2* 



i\\vor + \\zor) + [i + 



27^2 



b^ini 



(3.3) 



which is for later use. 

Let d = min{di,d3}. From ()2.25p and by Poincare inequality we can deduce that 



< 



d 
dt 

2|rfi-«2| e^rft||v(v(t) + 2(t))f+e^'^*(Ci(z;o,zo,i) + t^2(«o,2o,i)), ^ > 0. 



(3.4) 



di 
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Integrate (|3.4p to obtain 

min|l,/i ^) 

^ max{l,^}2M,-d,|^ re-^'^(*-)||V(z;(r)+z(r))f dr ^^'^^ 
«i Jo 

+ C'io(wo,-zo,t)! t > 0, 
where by (|2?22l) and ([Z!24|) we have 



Clo ("0 , ^0 , i ) = max{ l,fi} 



7^ \7c?i(i2 7^^ 
= max{l, /^} 

in which 



1 / M7«3 Jo fi^-^dd2d3 



f 8 ^ 2fc^ \ ^^^^^11^ ||2 ^ ||2) ^ / 46^ ^ 16q^ ^ Pb^ , 
\c?i7 fJ-jds J \'-f^ddid2 ^"^ddi ^j'^dd2d: 



•^0 [te-^'** < ^e-ie-^'^*/2^ if di = 2d2. 

Note that f3{t) — ;> 0, as t — > 00. On the other hand, multiplying p.l6p by e^''* and then integrating 
each term of the resulting inequality, we get 

I [ ^'"^^ {\H^W + UrW) dr + d2 e^'^^dl Vz;(r)||2 + || Vz(r)|n dr < ^^M^'"'- 

By integration by parts and using (|2.17p . we find that 

rf2^* e^'^WVvirW + ||Vz(r)f )dr < ^e'^'^' 



z-^'WHrn^ + \\z{t)r) iWM? + ll^of ) - / ide'^'^^iMrW + \\z{r)f)dT 



1 

~ 2 ^ 

< -^1 17| e^"^* + (||^;of + llzof ) + id f e^^'^^-^'^^^^Wv.r + \\z,f) dr + ffle'^'^* 
7" Jo 27^2 

< ^ + ^) h'me-"'' + (1 + 7^/3(06-^'') (Ikof + Ikof ), i > 0. 
Substituting p.7p into p.Sp we obtain that, for t > 0, 



(3.7) 
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max{ 1 , /i ^ } 



Wvitw + wmr < e--""' : ' {\\u^ + «o + + ^or + iiv^o + Mn + c^ioK, ^o, t) 

min|l,/i ^1 

2 max{l, — 



ax{l,/ 
min{l, 



--fdt 


'1 


G- 






.7 






+ I'D 


+ 


Wo 4 





max{l,A*} + — r /^WdkolP + ll^oin + 3 , + + ^ 



(3.8) 



Since (3{t) — and e as t 00, from p.8p we come up with 

limsup(||y(t)||2 + ||eW||2) <i?i, (3.9) 

where 

R -1 + max{l ^} + 16a2 ^ fc252 ^ 2|di -daP / 1 ^ ^\ |^| (3 

\73c?(iid2 72^6?! ^^^dd2d3 7^1 c?2 yd ^2 / / 

Moreover, for any given bounded set B C H there exists a finite time T^'^ > such that ||?/(i)|l2 + 
||^(t)||2 < i?i for any go £ B and ah t > T^'^. The combination of (|X2|) and (|3Jl) gives rise to 

hmsup(||w(t) + w(i)||2 + \\ip{t) + ^/j{t)\\^) 

|2 I I ././J.MI2 



= limsup \\y{t) - {v{t) + zm\' + \\ip{t) + mV) (3.11) 

< limsup (3(||y(0f + |k(t)f + ||z(t)f ) + UW) < 3{Ro+Ri). 

From the inequahty (12.37^ satisfied by p{t) = u{t) + v{t) — w{t) — z{t) and 9{t) = — i^it) and 
by Poincare inequality, similarly we get 

I (e^^*(|b(i)f + \\^^-'e(t)r) < ^^f^ e^''\Mv(t) - z{t))r ^3 ^2) 

+ e^-^* (Cg (wo , ^0 , i ) + (^'o , ^0 , i ) ) , i > 0. 

Integrate p. 121) to obtain 



min|l,/i ^} 

+ max{l, fi} (^^^ ^* e-^''^*'^) II V(«(r) - z(t)) f dr + Cn (z;o, ^o, 

where 



(3.13) 
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Using p.7p to treat the integral term in (I3.13p . we obtain 



\\p{t)f + \m)f < e'"'* rnn'^n (lh^o + - ^0 - zof + yo - M') 



< e 



V "1 

minjl, 



)(l|Vz;(r)f + ||V^(r)f ) dr + CnCt-o, ^o,i) 



(3.15) 



I-uq + Wo - wq - ^olP + Wfo - V'olP) + max{l, /i}Cii(wo, ^o, 



max{l, /^} 



2|di-rf2|' 

(ii(i2 











3^ 





62|r!| + (e-^'^*+7d/3(t))(||t;of + ||zo||' 



Again, since (3{t) and e-'^''* as t -> 0, from (jXTil and ([3J5)) we get 

|2 I l|/ir+M|2\ 



iimsup(ib(t)r + ii^wir) <^2, 



where 

i?2 = 1 + niax{l, /i} 



2|di - dap / I 
7c?id2 



\k + 2{Di-D2)\ 
Di 



(3.16) 

b^\n\. 

(3.17) 



Moreover, for any given bounded set B C H there exists a finite time T^' > such that ||]3(t)|p 
\\d{t)\\'^ < R2 for any go £ B and all t > T^'^ . The combination of ((X^ and gives rise to 



limsnp{\\u{t)-w{t)\\^ + y{t)~m\\') 

t—^OO 

= limsup {\\p{t) - (vit) - z{t))f + Mt) mf) 

<limsup [3(||p(t)f + \\vit)\\^ + ||z(t)f ) + ||0(t)||2] < 3(i?o + R2). 



(3.18) 



Finally, putting together ((3lT1) and (I3TT8)) and noticing ^OAO^i we reach the conclusion that 



limsup ||5(t)|P = limsup 



i II (u + w) + (?/ - w) f + 1 II (u + w) - (U - W) f 



(3.19) 



P + ll^-^f + ||«|p 



< limsup (||u + w|| + ||u - w|| +\\(p + ' 

< 7i?o + 3(i?i + i?2), 

and that for any given bounded set B C H there exists a finite time 

^B = max{T^'^T|■^TP''^}>0 

such that ||g(i)|p < Ki = IRq + 3(i?i + -R2) for any go e B and aU t > Tb- Therefore the lemma 
is proved with Ki = 7Ro + 3(i?i + -R2) in the description of an absorbing set Bo in (|3.ip and Ki is 
a constant independent of initial data. □ 

Next we show the absorbing properties of the {v, z) components of the extended Brusselator 
semiflow in the product Banach spaces [L^P(ri)]^, for any integer 1 < p < 3. 
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Lemma 3.2. For any given integer 1 < p < 3, there exists a positive constant Kp such that the 
absorbing inequality 

limsup \\iv{t),zit))\\%, < Kp (3.20) 

is satisfied by the {v,z) components of the extended Brusselator semiflow {S(t)}t>o for any initial 
data go G H. Moreover, for any bounded set B C H there is a finite time > such that 
\\{v{t), z{t))\\^j^2p < Kp for any g^ E B and all t > Tg. 

Proof. The case p = 1 has been shown in Lemma 13.11 According to the solution property (|2.ip 
with Tmax = oo for all solutions, for any given initial status go E H there exists a time to £ (0, 1) 
such that 

S{to)gQ eE = [H^{n)]^ ^ LS(f^) L''(f^). (3.21) 
Then the weak solution g(t) — S{t)go becomes a strong solution on [to, oo) and satisfies 

Si-)go e C([io, oo); S) n L^{to, ^; n) C C{[to, ^); L'^(r!)) C C([io, (x,);L\n)), (3.22) 

for n < 3. Based on this observation, without loss of generarity, we can simply assume that 
go e L^(r2) for the purpose of studying the long-time dynamics. Thus parabolic regularity p.22|) 
of strong solutions implies the S{t)go E E C L^(f7),t > 0. Then by the bootstrap argument, again 
without loss of generality, one can assume that go G n C L^(ri) so that S{t)go S 11 C L^(r2),< > 0. 
Take the inner-product ( (|1.2p . v^) and ( (jl.Sp . z^) and sum them up to obtain 

~ (IkWIlia + ll^(i)llie) + 5rf2 {h{t?^v{tW + \\z{t?Vz{t)r) 

(bu{t,x)v^{t,x) - u^{t,x)v^(t,x) + bw{t,x)z^{t,x) - w^it,x)z^{t,x)) dx (3.23) 



+ D2 [{z{t,x) - v{t,x))v^{t,x) + {v{t,x) - z{t,x))z^{t,x)] dx. 
Jn 

By Young's inequality, we have 

[{buv^ - u^v^) + {bwz^ - w^z^)] dx<^(^J b^{v* + z^) dx- J {u^v^ + w^z^) dx 



and 



[(z — v)v^ + (v — z)z'^'^ dx < 



6 6 / V6 6 



dx = 0. 



Substitute the above two inequalities into (I3.23P and use Poincare inequality, we get the following 
inequality relating ||(w,z)|||6 to ||(w,z)||^4, 

I iMrnh + umh) + ^7^2 {\\vmie + imwie) 

* fn (3-24) 

< ^ {Hmh + \\zit)\\%) + -d2 iwvv'm' + iivz3(i)f ) < 36^(ik(t)iii. + \\zmu). 

Similarly by taking the L inner- product (([112]), u^) and ((IT3t.z^) we can get the corresponding 
inequality relating ||(w,z)||^4 to ||(w,z)|p, 

^ {Mt)\\u + \\zmu) + 37d2 (ii«(t)iii4 + \\zit)\\i,) 

(3.25) 

< I {\Ht)\\h + wmih) + 3d2 {\\yvHt)r + \\\7zHt)r) < 2bH\\v{t)r + 
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Applying Gronwall inequality to inequalities p.25p . p.24p . and using (|2.17p . we get 

\Ht)\\l. + \\z{t)\\l. 

<e 







<e-^'^^*Ci2 (||«o|li« + IkolliO + t > 0, 

where C12 is a positive constant, and then 

\\vmie + \\zm% 







<e 



<e-'^'^=^*Ci3 (Ikollie + llzollie) + t > 0, 

where C13 is a positive constant. It follows that 

hmsup (Ib(t)lli. + ||z(t)||i.) <K, = l + (3.26) 

hmsup (||^;(^)||ia + ||z(t)||ie) < i^s = 1 + (3-27) 

Thus p.20p is proved and the last statement of this lemma is also proved. □ 



4. Asymptotic Compactness 

In this section, we show that the extended Brusselator semiflow {S(t)}t>o is asymptotically 
compact through the following two lemmas. Since H^ip.) M> L'^{U,) and Hq{^) ^ L^{il) are 
continuous embeddings, there are constants (5 > and 77 > such that || • \Wi < <5||V(-)|p and 
II • Hie < '7||V(-)P- We shall use the notation ||(yi, j/2)|P = WviW + Wv^W^ for conciseness. 

Lemma 4.1. For any given initial data go £ Bq, the {u,w) components of the solution trajectories 
g{t) ^ 5(t)5o of the IVP (lllil) satisfy 

||V(M(t),u;(t))f < Qi, fort>ri, (4.1) 

where Qi > is a constant depending only on Ki and hut independent of initial data, and 
Ti > is finite and only depends on the absorbing ball Bo. 



Proof. Take the inner-products ( (jl.ip . —Au{t)) and ( (|1.4p . —Aw{t)) and then sum up the two equal- 
ities to obtain 
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~\\viu,w)r+d,\\Aiu,w)r+ib+k)\\viu,w)r 



= a{Au + Aw) dx — / {u'^vAu + w'^zAw)dx — N / {ipAu + ipAw) dx 
Jn Jn Jn 

- Di {\Vu\^ -2Vu-Vw +\Vw\^)dx 
Jn 



^ ■ 4 ' 4 ' 2 
It follows that 



-\\Viu,w)f + 2ib + k)\\V{u,w)f 

< |. ^a'\n\ + iv2||(^,^)f ) + 1 iWuYWvr + ll-^lPll-f ) 

< |- {a^\n\ + iv^iK^, ) + J (lkinivz.||4 + ||z|P||V^«||4) , t>0. 



(4.2) 



By the absorbing property shown in Lemma |3.1[ there is a finite time Tq = Tq{Bq) > such that 
S{t)Bo C Bo for all t > Tq. Therefore, for any go G Bq, by ((3l1) we have 



ll(«W,^"W)f + ll(«W,^W)f + < i^i, for t > To. 

Substitute (l43l) into (1421) to obtain 



|||VKz«)f <|||VKz«)f +2(fe + fc)||VK^)f 



< 



\\yiu,w)r 



di 



'■\n\+N^Ki), t>To, 



which can be written as the inequality 



±<ap+^{a'\n\+N'K,), t > To, 
at di 



where 



p{t) = \\\7{u{t),w{t))\\'^ and a{t) 



di 



Pit). 



In view of the inequality (|2.20p , (I3.3p and (|4.3p , we have 



t+i 



\\Vy{T)\\^dT< 



2\di-d2\ 



2 ft+l 



d\ 



2k 



\\V{v + z)rdr+- I ||S(r)fdr 



t+i 



d[ ( "1^(^)11^ 



-{Mt)^ + \\z{T)f + 2a'M)dT] <Ci4, for t > To, 
t 1 



where 



Ci4 — 



4|di -dzl 
dld2 



27d2 



+ — ( (4fc^-i + + -{Ki + 2a^\n\) 



(4.3) 



(4.4) 



(4.5) 



(4.6) 
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From the inequality (|2.32p . p.3p and (|4.3I) and with a similar estimation, there exists a constant 
Ci5 > such that 

pt+i 

1 1 Vp(t) f dr < Ci5 , for t > n. (4.7) 

According to (|2^ . we can put together ((33)) . (|4?6t and (|4J| to get 

r-t+l /-t+l 

p(T)dr= / (||Vu(r)f + ||V«;(r)f)dT 



(4.8) 



1 

< - 
- 2 



(l|V(2/(r) - (vir) + z(r))f + ||V(p(r) - («(r) - z(r))f ) dr 

t+i 

2 I I|V7„/_M|2 I ||V7,'„, I ,M|2 I ||v7/'„, .,M|2^ 



< / (l|V2;(r)f + |lVp(r)f + ||V(« + z)f + ||V(« - z)f ) dr 



4 

< Ci4 + Ci5 + — 
^2 



='Ci6, forOTo. 



Now we can apply the uniform Gronwall inequality, cf. [341 Lemma D.3], to (j4.5p and use (|4.8p to 
reach the conclusion (14. ip with 

gi = (^Ci6 + |- {a'\n\ + N^K,)^ eS^K,C,e/d, 

and Ti = To(So) + 1. The proof is completed. □ 

Lemma 4.2. For ant/ given initial data go £ Bq, the {v,z) components and (ifjip) components of 
the solution trajectories g{t) — S{t)gQ of the IVP (jl.l4l) satisfy 

||V(«(t),z(t))f + ||V(^(t),V(t))f <Q2, fori>T2, (4.9) 

where Q2 > is a constants depending on Ki and but independent of initial data, and T2 (> 
Ti > 0) is finite and only depends on the absorbing ball Bq . 

Proof. Take the inner-products ( (|1.2p . —Av{t)) and ( (jl.51) . — Az(t)) and then sum up the two equal- 
ities to obtain 

~\\Viv,z)f + d2\\A{v,z)\\^ = - [ b{uAv + wAz)dx 

+ / {u'^vAv + w'^zAz) dx- D2 / [{z - v)Av + {v ~ z)Az] dx 
Jn Jn 

<^\\Aiv,z)\\^ + ^\\iu,w)\\^ + [ {u^"" + w^z'')dx (4.10) 
^ "2 «2 Jn 



D2 f {\Vvf -2\/v\/z+\Vzf)dx 
Jn 



<^||A(«,z)f + .^||(«,u,)f + 1 / {u''v^ + w''z^)dx, t>0. 
2 "2 «2 Jn 



Since 



\\viv{t),z{t})r = ^{{v,Av) + {z,Az)) < i {\\v{t)r + \\z(t)r + w^vitw + 
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in (|4.10p we have 

^\\Aiv,z)r >d,\\viv,z)r -^\\iv,z)r. 

Then by using Holder inequahty and the embedding inequahties mentioned in the beginning of this 
section and by Lemma 14711 from the above inequahty (|4.10|) we get 

|||V(z;,z)f <|||VKz)f + 2d2||VKz)f 

<d2\\{v,z)r + ^Uu.wW + ^(Ikllialkllia + \\w\\U4h) 



- ("^^ + ¥) + + \\Wwf)\\V{v, z) 



\\y{v,z) 



t > Ti. 



Applying the uniform Gronwall inequality to 
dt 



||V(t;,z)f <^||V(^,z)f + 



262 
d2 



t>Ti, 



and using p.3p . we can assert that 

||V(«(t),^(t))f <Ci7, for t>T, + l, 

where 



(4.11) 
(4.12) 



1 



1 



2jd2 



b^\n\ 



2b' 



^2rj''Ql/d2 



Next take the inner-products ( (|1.3p . ip(t)) and ( (|1.6p . tbd)) and then sum up the two equalities 
to get 



i|||(^»f + d3|iV(^»f + (A + N)\\i^,i;)r < ^ {WMf + ||(^, V)f ) 



so that 



t+i 



||V(^(T),V'(r))f dr < — I|(^(i),V(0)f 



k 

'2d3Jt 



(||(u(r),u;(r))f + ||(^(r),V'(r))f)dr< 



(1 + k)Ki 
2d. 



(4.13) 



t>To. 



Then take the inner-products ( ^.31 . —Au^d)) and ( (11.61) . —Aibit)) and sum up the two equalities 
to obtain 

^ ||V(^, + ds\\A{^, + (A + N)\\S/{^, 



2dr 



k{uA(p + wAijj) dx ~ D3 / [(V- - v)Aip + [ip - -4')A^p] dx 



<d4A{^,^)f + —{\\u\\'' + \\wf)~D^ / (|V^|2-2V^- W^+IW^ndx 



4c(3 
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SO that 

-^||V(^»f + 2(A + 7V)||V(^,^)f <^i^i, t>To. (4.14) 
at 20,3 

We can apply the uniform Gronwah inequahty to (|4.14p with the aid of (|4.13p to reach the estimate 
||VMt),V(i))f <^(l + fc + fc')e'(^+^) ='Ci8, t>To + l. (4.15) 
Finally let Q2 = Cn + Cig and T2 = Ti + 1. By (|iT^ and dHH]), we see that (gH) is proved. □ 

5. Global Attractor and Its Properties 

In this section we reach the proof of Theorem 11.51 (Main Theorem) on the existence of a global 
attractor, which will be denoted by £/, for the extended Brusselator semiflow {S{t)}t>o and we 
shall show several properties of this global attractor £/, namely, the regularity of £/, the property 
of being an {H, E) global attractor, and the finite Hausdorff and fractal dimensionality. 

Proof 0/ Theorem 1. In Lemma [3.1l we have shown that the extended Brusselator semiflow {S{t)}t>o 
has an absorbing set Bq in H. Combining Lemma |4. II and Lemma 14.21 we proved that 

\\S{t)go\\l <Qi+ Q2, for t > T2 and for go e Bq, 

which implies that {S{t)BQ : t > T2} is a bounded set in E and consequently a precompact set in H. 
Therefore, the extended Brusselator semiflow {S{t)}t>o is asymptotically compact in H. Finally 
we apply Proposition 11.61 to reach the conclusion that there exists a global attractor £/ in H for 
this extended Brusselator semiflow {S{t)}t>o. □ 

Now we show that the global attractor jz/ of the extended Brusselator semiflow is an {H, E) 
global attractor with the regularity C L°°(i7). The concept of {H,E) global attractor was 
introduced in [3]. 

Definition 5.1. Let {S(t)}f>o be a semiflow on a Banach space X and let F be a compactly 
imbedded subspace oi X. A subset yi of F is called an {X, Y) global attractor for this semiflow if 
A has the following properties, 

(i) .A is a nonempty, compact, and invariant set in Y . 

(ii) A attracts any bounded set B C X with respect to the K-norm, namely, there is a time 
T = t{B) such that E(i)B C F for t > t and distY{^{t)B,A) 0, as t ^ 00. 

Lemma 5.2. Let {gm\ be a sequence in E such that {gm} converges to go d E weakly in E and 
{grn} converges to go strongly in H , as m ^ 00. Then 

lim S{t)g„i ~ S{t)go strongly in E, 

and the convergence is uniform with respect to t in any given compact interval [toi^i] C (0,oo). 

The proof of this lemma is seen in j45l Lemma 10]. 

Theorem 5.3. The global attractor £/ in H for the extended Brusselator semiflow {S{t)}t>o is 
indeed an {H,E) global attractor and si/ is a bounded subset in L°°(f2). 

Proof. By Lemma 14.11 and Lemma 14.21 we can assert that there exists a bounded absorbing set 
Bi C E for the extended Brusselator semiflow {S{t)}t>o on H and this absorbing is in the i?-norm. 
Indeed, 

Bi = {geE:\\g\\l = \\\7g\\^<Qi+Q2}. 
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Now we show that this scniiflow {S{t)}t>Q is asymptotically compact with respect to the strong 
topology in E. For any time sequence {in},tn — > oo, and any bounded sequence {gn} C E, there 
exists a finite time ta > such that S(t){gn} C Bq, for any t > Iq. Then for an arbitrarily given 
T > to +T2, where T2 is the time specified in Lemma 14. 2| there is an integer hq > I such that 
tn > 2T for all n > uq. According to Lemma [4. II and Lemma [4.21 



{S{tn — T)gn}n>no is a bouudcd set in E. 

Since E is a Hilbert space, there is an increasing sequence of integers {nj}°°^i with ni > no, such 
that 

lim S{tn. - T)g„- = g* weakly in E. 

By the compact imbedding E ^ H, there is a further sebsequence of {nj}, but relabeled as the 
same as {%■}, such that 

lim S{tn^ - T)gn, g* strongly in H. 
Then by Lemma I5.2[ we have the following convergence with respect to the E'-norm, 
lim 5(t„.)5„. = lim S{T)S{tn^ - T)gn, = S[T)g* strongly inE. 

This proves that {S{t)}t>o is asymptotically compact in E. 

Therefore, by Proposition 11.61 there exists a global attractor s^e for the extended Brusselator 
semiflow {S{t)}t>o in E. According to Definition 15.11 and the fact that Bi attracts Bq in the E- 
norm due to the combination of Lemma [4. II and Lemma [4.21 we see that this global attractor s/e 
is an [H, E) global attractor. Moreover, the invariance and the boundedness oi H and in E 

imply that 

£^E attracts ^ in E, so that si/ C £/e\ 
s/ attracts s/e in i?, so that s/e C 

Therefore, s^/ = s/e and we proved that the global attractor jz/ in _ff is indeed an E) global 
attractor for this extended Brusselator semifiow. 

Next we show that jz/ is a bounded subset in L°°(ri). By the {Lf ^ L°°) regularity of the analytic 
Co-semigroup {e'^*}t>o stated in [Ml Theorem 38.10], one has e^* : LP(rj) — > L°°(f7) for t > 0, 
and there is a constant C{p) > such that 



e 



^*\\L{hp,i.'=-) < C{p)t ^, t>0, where 71 = dim 17. (5.1) 

By the variation-of-constant formula satisfied by the mild solutions (of course strong solutions) , for 
any g G jz/ (c E), we have 



\\S{t)g\\L^ < ||e^*|U(,..,i.)||5||+ / \\e^^'-''^hiL^L^)\\f{S{a)g)\\da 

,t " (5-2) 

<C(2)t-t||5||+ / C{2){t-ariL{Q,,Q2)\\S{a)g\\Eda, t>0, 
Jo 

where C(2) is in the sense of (|5.ip . and L{Qi,Q2) is the Lipschitz constant of the nonlinear map 
/ restricted on the closed, bounded ball centered at the origin with radius Qi + Q2 in E. By the 
invariance of the global attractor surely we have 

{S{t).s/ : i > 0} C So (C i/) and {Sit)^/ : t > 0} C Bi (c E). 
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Then from we get 

\\S{t)g\\L^ <C{2)Kit-i + f C{2)L{Qi,Q2){Qi + Q2){t~<j)-Ua 

Jo (5.3) 

= C{2)[Kit-^ + AL{Qi,Q2){Qi + Q2)ti], ioi t > 0. 
Specifically one can take i = 1 in (|5.3p and use the invariance of to obtain 

Ml^ < C{2){Ki + 4L(Qi, g2)(Qi + Q2)), for any <? e 
Thus the global attractor ^ is a bounded subset in L°°(f7). □ 

Next consider the Hausdorff and fractal dimensions of the global attractor . The background 
concepts and results can be seen in [36l Chapter V]. Let q,n — \iuYsn\)^^^ qm{t), where 

q™(t)=sup sup / Tr[(^ + /'(5(r)5o))r„(T)] drV (5.4) 

go&^ gieH,\\gi\\=i V r Jo / 

i—l,--- ,m 

in which Tr [(yl + /'(S'(T)go))r™(T)] is the trace of the linear operator + /'(S'(T).go))r,„(T), f'{g) 
is the Frechet derivative of the Nemytskii map / in (|1.14p . and Tm{t) stands for the orthogonal 
projection of the space H on the subspace spanned hy Gi{t), ■ ■ ■ , Gm{t), with 

G.(t) =i(5(0,5o)g<, i^l,---,m. (5.5) 

Here f'{S{T)go) is the Frechet derivative of the map / defined by (|1.13p at 5(t)(7o, and L{S{t),gQ) 
is the Frechet derivative of the map S{t) at g^, with t fixed. 

The following proposition, cf. [JB^, Chapter 5], will be used to show the finite upper bounds of 
the Hausdorff and fractal dimensions of this global attractor s/. 



Proposition 5.4. // there is an integer m such that g.,„ < 0, then the Hausdorff dimension dni-B/) 
and the fractal dimension dpi-s^) of satisfy 

dni^) < m, and dpi^^) < m max (1 + ) < 2m. (5.6) 

l<j<m-l \ \q,n\ J 

It can be shown that for any given t > 0, S{t) is Frechet differentiable in H and uniformly 
Frechet differentiable in Mf. Its Frechet derivative at go is the bounded linear operator L{S{t),gQ) 
given by 

L{S{t), go)Go = G{t) = {U{t), V{t), m,W{t), Z{t), ^{t)), 

for Go = (C/o, Vo, $0, Wo, Zo, ^0) G where {U{t),V{t),<^>{t), W{t), Z{t),^{t)) is the weak solution 
of the following extended Brusselator variational equation 
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dlJ 

— = diAU + 2u{t)v{t)U + u^{t)V - (6 + k)U + Dj_iW - U), 
dV 

— = rfaAy - 2u{t)v{t)U - u^{t)V + bU + Dgf^ - V), 
at 

— = d3A$ + fcC/ ~ (A + N)^ + D3(* - $), 

^ = diAM^ + 2w(i)z(t)H/ + w^{t)Z -{b + k)W + Di{U - W), (^■'') 

— = rfaAZ - 2w{t)z(t)W - w^(i)Z + 6W^ + DafV - 
at 

9* 

— = dgA* + fcVK - (A + A^)* + L>3(<I> - 
at 

{/(O) = C/o, 1/(0) = Fo, $(0) = $0, VF(0) = Wo, Z(0) = *(0) = *o- 

Here g(t) ~ {u{t),v{t), ip{t),w{t), z{t),ilj{t)) — S{t)go is the weak solution of (|1.14p with the initial 
condition ^(0) = go. The initial value problem (j5.7l) can be written as 

rlC 

— = {A + f{S{t)go))G, G(0) = Go. (5.8) 
As we have shown, the invariance of £/ implies £/ C Bq f) Bi, so that 

sup \\S{t)go\\^ < Ki and sup \\S{t)go\\l < Qi + Q2- (5.9) 

Theorem 5.5. The global attractors £/ for the extended Brusselator semiflow {S{t)}t>o has a 
finite HausdorfJ dimension and a finite fractal dimension. 

Proof. By Proposition [53J we shall estimate Tr [{A + f {S{T)go))Tm{T)]. At any given time r > 0, 
let {Cj{T) : j = 1, • • • , m} be an i/-orthonormal basis for the subspace 

r„,(r)ff = Span{Gi(r),... ,G,(r)}, 

where Gi(i), • • • , Gm{t) satisfy (|5.8p and, without loss of generality, assuming that the initial vectors 
Gi^o, ■ ■ ■ , Gm^o E^re linearly independent in H. By the Gram-Schmidt orthogonalization scheme, 
CjiT) G E and Cji''') is strongly measurable in T,j = 1, ■ • ■ Let c?o = min{c?i, ^2, ds}. Denote 
the components of Qir) by Cj(T), i = 1, • ■ • , 6. Then we have 



Tr [{A + /'(5(r)5o)r„(r)] = ^ ((AO(r), 0(r)) + (/'(5(r)5o)0(r), 0(r))) 

m 

< -doJ2 ||VO(T)f + J1+J2 + J3, 



where 



Ji^J2 MrHr){\(:j{T)\^-Cj{T)(:]{r))dx 



+ J2 / 2z.(r)z(r)(|C|(r)|2-Cl(r)C|(r))rfa:, 
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<E / ("'(T)lCl(r)||C^(r)|+«;2(r)|C|(r)||C|(r)|)dx, 



and 



•^3 = E / (-(& + fc)(ICl(r)P + |C^(r)p)+fe(Cj(r)C|(r) + C^(r)C|(r)))dx 



+ E / (fc(Cj(r)Cj'(r) + Cl(r)Cj*(r))-(A + iV)(|C|(r)p + |C^(r)n)dx 

.-1 -'O 



E / (^i(C}W-C»)%i?2(C|(r)-C|(r))%i^3(Ci'(r)-Cf(r))^) 



^E( lj{Q]{r)C^{r)+e,{T)C^{T))dx+ lj{Cj{r)Cf{T)+Cf{r)Cj{T))dx 



n 



By the generalized Holder inequality, the embedding i?o(r2) ^ L'*(r2) (for n < 3) and (|5.9p . we get 

m 

Ji<2Ell '"('r)IU''lk('r)IU* (IICj (■r)llL-i + IIC}(t)IIl*IICj 



-2E |k(r)|U.||z(r)|U4 (llC^(r)lli. + ||C^(r)|U4||Cj^(r)|U.) 

(5.11) 



<4E l|5(r)golli4||0(r)||i. < 45E II V5(r)goiniO(r) 



li. 



<4(5(Qi + g2)EllO(T)|li4, 

where 5 is the embedding coefBcient given in the beginning of Section 4. Now we use the Garliardo- 
Nirenberg interpolation inequality for Sobolev spaces [Ml Theorem B.3], 

IICIIw^,. < CIICII^.,, llClli"', for C e W^'^in), (5.12) 
provided that p, q,r > 1,0 < 9 < 1, and 

Ti I Tl\ Ti 

k <9\m — (1 — 0) — , where n — dinifi. 

p \ q J r 

Here with W'''P{fl) = L^{n),W"^-'^in) = H^{n), L'' (ft) = L^ifl), and 61 = n/4 < 3/4, it follows 
from ([51^ that 

|IO(r)|U4 <q|VO(T)||t||C^.(r)||i-t =q|VO(r)||t, j- = l,... ,m, (5.13) 
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since ||Cj(''')ll — Ij where C is a universal constant. Substitute (|5.13p into (|5.1ip to obtain 

m 

Ji <4<5(Qi + Q2)C2^ I|VO(t)||*. (5.14) 

Similarly, we can get 

m m 

J2 <<5(Oi+Q2)5]|IO(r)||i4 <5(Qi + Q2)C2^ l|VO(T)||t. (5.15) 
Moreover, we have 

m 

J3<^(fe + fc)||0(r)f =m(6 + fc). (5.16) 
Substituting ((gUl) . and (j^TTBl) into we obtain 



Tr[(A + /'(5(r).9o)r„,(r)] < -do^||VO(T)| 



(5.17) 



+ 5(5(gi + Qs)^' I|V0(t)||5 + m(6 + fc). 

By Young's inequality, for n < 3, we have 

5<S(Qi + Q2)C2 E II^OWII^ < y E l|VO(r)f + mQ3(n), 

where Q3(?i) is a universal positive constant depending only on n = dim (J7). Hence, 

, m 

Tr [(A + /'(5(r).9o)r,„(r)] < ^ II^O(T)f + Tn{Q^{n) + 6 + fc), r > 0, go e ^. 

According to the generalized Sobolev-Lieb-Thirring inequality Appendix, Corollary 4.1], since 
{Ci(''')r'' jCm(''')} is an orthonormal set in there exists a universal constant Q* > only 
depending on the shape and dimension of $7, such that 



Therefore, 

TrP + /'(5(T)go)r™(T)] < -^^mi+^ +m(Q3(ri) + 6 + fc), r > 0, go e ^. (5.18) 

Then we can conclude that 

= limsup q^it) 

= limsup sup sup (I f Ti-[{A + f{S{T)go))T„^(T)]dT] 

t^oc gaesi g,e-ff,||g, 11 = 1 K'^Jo J (5.19) 

z=l,--- ,m 

< - ^^mi+* + m{Qz{n) + b + k)<0, 
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if the integer m satisfies the following condition, 

^-l<^J^M_i_ \n\<m. (5.20) 

According to Proposition l5.4[ we have proved that the HausdorfF dimension and the fractal dimen- 
sion of the global attractor £/ are finite and their upper bounds are given by 

dni'S^) < m and c?i?(^) < 2to, 

where m satisfies (j5.20p . □ 

As a conclusion remark, in this work it is shown that global dynamics of this extended Brusselator 
system of cubic-autocatalytic, partially reversible reaction-diffusion equations with linear coupling 
between two compartments and with homogeneous Dirichet boundary condition is dissipative and 
asymptotically determined by a finite-dimensional global attractor in the Hq product phase space. 

We emphasize that even though this multi-dimensional reaction-diffusion system lacks the math- 
ematical dissipative condition, the existence of a global attractor is established without assuming 
any non-negativity of initial data (or solutions) and without any conditions imposed on any of 
the positive parameters. The contributed methodology of grouping and re-scaling estimation suc- 
ceeded in proving the dissipative longtime dynamics has a promising potentiality to be adopted and 
adapted for future investigations of many open problems and complex coupling models in system 
biology and in network dynamics. 

The extensions of the obtained results can be pursued in several directions. Semicontinuity of the 
global attractor with respect to one of the diffusive coefficients and/or with respect to the reverse 
reaction rate coefficient when it tends to zero is a meaningful question of singular perturbation of 
global dynamics. One can also consider similarly coupled reaction-diffusion systems on a higher- 
dimensional domain of space dimension n > 3 and on an unbounded domain to work with various 
different phase spaces. 

A challenging problem is to study the two-cell or multi-cell coupling system of the Brussela- 
tor equations or the Gray-Scott equations in which one subsystem is defined on a subdomain and 
the other subsystem is defined on another non-overlapping subdomain with the flux-type coupling 
through the joint boundary, such as in many biological and physiological models of signal transduc- 
tions. This kind of reaction-diffusion systems of investigation will be very close to the real-world 
problems in cell biology and in physical chemistry. 
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